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Abstract. We make the case that the Coulomb- plus linear quark confinement potential predicted by 
lattice QCD is an approximation to the exactly solvable trigonometric Rosen- Morse potential that has the 
property to interpolate between the Coulomb- and the infinite wells. We test the predictive power of this 
potential in the description of the nucleon (considered as a quark-diquark system) and provide analytic 
expressions for its mass spectrum and the proton electric form factor. We compare the results obtained 
in this fashion to data and find quite good agreement. We obtain an effective gluon propagator in closed 
form as the Fourier transform of the potential under investigation. 

PACS. PACS-key 19.39.Jh (non-relativistic quark models), 13.40.Gp (electromagnetic form factors) 



1 The quark potential from lattice QCD. 

The strong interactions of quarks, the fundamental con- 
stituents of hadrons, are governed by the Quantum Chro- 
modynamics (QCD), the non-Abelian gauge theory with 
the gluons as gauge bosons. As a consequence of the non- 
Abelian character of QCD, the quark interactions run 
from one- to many gluon exchanges over gluon self-int- 
eractions, the latter being responsible for the so-called 
quark confinement, where highly energetic quarks remain 
trapped but behave as (asymptotically) free particles at 
high energies and momenta. The QCD equations are non- 
linear and complicated due to the gluon self-interaction 
processes and their solution requires employment of highly 
sophisticated techniques such as discretization of space 
time, so-called lattice QCD. Lattice QCD calculations have 
established themselves as a reliable tool for non-perturba- 
tive analysis of QCD. The outcome (in the quenched ap- 
proximation) is a linear confinement potential with energy 
increase, be it quark-anti-quark (QQ), 
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(1) 



or two-body potential between quarks (so called A type) , 



V3Q(ri,r 2 ,r 3 ) = -A A ^2 



r t - r,- 



Kj 



v 3 \+C A , 
. (2) 

in obvious notations. The Coulomb-like piece is associated 
with short range one-gluon exchange in the perturbative 
regime, whereas the long-range linear part relates to non- 
perturbative effects and is attributed to flux-tube QQ, or 
QQ links. Its strength is then associated with the respec- 
tive string tension [I] . Detailed analysis of the values of the 



constants of Vqq and V^q has been performed in Ref . [2] . 
The three-quark (3Q) Coulomb- plus linear potential, or, 
versions of it, has found repeatedly application to baryon 
spectroscopy 3,4 . In view of this, its generalization to an 
exactly solvable potential is of interest. 

We begin with first drawing attention to the proxim- 
ity of the two-body Coulomb- plus linear potential (be it 
for QQ, QQ, or Q(QQ) systems), to (— cot|r|). Indeed, 
this is immediately seen from the corresponding Taylor 
expansion, 

_cot|r|«-ir + i|r|. (3) 

This expression shows that the absolute values of the stre- 
ngths of the linear to Coulomb potentials are in ratio 1:3, 
a value that fits quite reasonably into the range of the 
a A '■ A a ratio of approximately 1/4 — 1/2 reported by the 
lattice QCD analysis [2]. In fact, (— cot |r|) is part of the 
more general and exactly solvable trigonometric potential 



VtflAf (|z|) = — 2&cot |z + a(a + 1) esc 2 |z| 



d' 



(4) 



known in supersymmetric quantum mechanics (SUSYQM) 
under the name of the trigonometric Rosen-Morse poten- 
tial [5] and displayed in Fig. 1. The potential in eq. (0} 
is given in terms of dimensionless argument and parame- 
ters. The figure shows that v t RM interpolates between the 
Coulomb- like potential (associated with the one-gluon ex- 
change) and the infinite well. The latter provides the ad- 
equate scenario for the asymptotic freedom in so far as it 
describes trapped but free particles. The intermediate re- 
gion of the by and large linear confinement potential can 
be attributed, as usual, to gluon flux tube links between 
the quarks. The length scale, d, will be considered as a 
free parameter to fit data. The Taylor expansion of the 
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Fig. 1. The trigonometric Rosen-Morse potential (solid 
line) and its proximity to the Coulomb- plus linear potential 
(dashed line) for the toy values a = 1, b = 50 of the parameters. 



esc 2 term reads, 



a(a + 1) esc 2 |z| 



a(a + 1) a(a + 1) 
z 2 + 15 1 



(5) 



The parameter a can be viewed as relative angular mo- 
mentum in the case when it takes integer non-negative 
values, a parametrization used in ref. [BJ. If so, the esc 2 
term acquires meaning of a non-standard centrifugal bar- 
rier. 

We were able to trace back VtRM to ref. [7] where it 
apparently has been found for the first time as a pure 
mathematical construct while completing the SUSYQM 
compilation of shape invariant potentials that are exactly 
solvable by means of the hypergeometric differential equa- 
tion. The VtRM potential is therefore one of the new poten- 
tials predicted by SUSYQM [5]. We here suggest that the 
QQ interactions are ruled by the following trigonometric 
quark confinement (TQC) potential 

K,q QC (zi,z 2 ,z 3 ) = -26^cot|z 4 -zj| 



+a(a + 1) ^ esc 2 |z 4 - Zj|, z t 



r, 



(6) 



The great advantage of the TQC potential over the Cou- 
lomb- plus linear potential (compared in Fig. {T])) is that 
while the latter is neither especially symmetric, nor ex- 
actly solvable, the former is both, it is exactly solvable, and 
has in addition the dynamical 0(4) symmetry [BJ for a = 
I = 0, 1, 2, and 50(2, 1) for more general a values [8]. 



2 The exact single particle basis of v tRM - 

In this section we briefly review for the sake of complete- 
ness of the presentation the exact solutions of the three di- 
mensional single particle Schrodinger equation with VtRM-, 



as it appears in the quark-diquark approximation to nu- 
cleoli structure [BJ[9] . From here onward we shall identify 
the a parameter with the relative quark-diquark angular 
momentum. The above equation is solved in the usual way 
in separating variables using the ansatz 



K(Ki)(\z\) 



(7) 



The exact solutions of the radial part have been con- 
structed in ref. [BJ on the basis of the one-dimensional 
solutions found in [5] as: 



K 



(Kl)[ 



= A^sin^lzle-^W* 2 * 1 



# = 0,1,2,. 



{K+l)) , , 
K-l (C0t|z|), 

1 = 0,1, ...,K. (8) 



Here, Nfj^l) i s a normalization constant and z is the rel- 
ative quark-diquark distance. The i?£ a '^(cot |z|) functions 
are the non-classical Romanovski polynomials [HHHT] which 
are defined by the following Rodrigues formula, 



-13+1 



d" 



— e" QCOt " x (l + x 2 )P- 1+n , (9) 



dx 



where x — cot |z| (see ref. [12; for a recent review). The 
energy spectrum of the TQC potential is calculated as 



1 



(K+l) 



+ (K + 1) 2 , K = n + l, (10) 



where n = 0,1,2,.... Therefore, particular quark levels 
bound within different TQC potentials (distinct by the 
values of I) happen to carry same energies and align to 
series of states (multiplets) characterized by the supe- 
rior quantum number K (see Fig. 2). One easily recog- 
nizes that the energy in eq. (TTTTj) is defined by the Balmcr 
term and its inverse of opposite sign, thus revealing 0(4) 
as dynamical symmetry of the problem. Correspondingly, 
the i^-series belong to the irreducible 0(4) representa- 
tions of the type (y > t) an( ^ ^ nas ^ ne meamn g of four- 
dimensional angular momentum. Upon coupling the quark 
spin to them, one finds the reducible 0(4) multiplct ( y> y) 

(|, 0) (0, ^), which consists of the K parity mates \ , 
■■■,(K — , and the one unpaired state of maximal spin 
J max = K + i of either positive, or, negative parity. 



2.1 The nucleon spectrum. 

Comparison of the TQC spectrum in Fig. 2 to the nu- 
cleon excitations reported by [13] reveals an amazing coin- 
cidence. All the observed nucleon resonances with masses 
below 2.5 GeV do indeed fall into K — 1,3,5 multiplets 
from which only the Fn and H i t n states are still "miss- 
ing", an observation due to refs. [14] . The K = 2, 4 lev- 
els have been attributed to entirely "missing" resonances. 
The unnatural parity of the K = 3, 5 states would re- 
quire to account of the internal structure of the diquark. 
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Fig. 2. Schematic presentation of the lowest K levels in eq. (JHJ. 
They unite states of same energies from different TQC poten- 
tials, here denoted by vi, whose angular momenta vary accord- 
ing to the rule I = 0,1, 2, K. 



This phenomenon can be interpreted as dominance of a 
quark-diquark configuration in nucleon structure (see last 
reference in [TJ] for more details) meaning that (i) eq. (10) 
can be directly employed to fit data on the nucleon spec- 
tra, (ii) matrix elements of transition operators can be 
evaluated in the basis of eq. (8). The nucleon spectrum is 
fitted by the following potential parameters [6] 



6 = 5.85, d=2.31fm /j, = 1.06 fm~ 



(11) 



where // stands for the reduced mass of the quark -diquark 
system. The wave function of the nucleon ground state, 
K = 0. 1 = 0. is obtained as 



Tl m (\z\) =N e~ b ^ sin 



No 



4b(b 2 



1) 



1 



= -2tt6 



(12) 



In order to illustrate the predictive power of the TQC po- 
tential and the efficiency of its exact single particle basis 
we shall exploit in the next section the potential param- 
eters in eq. pT|) and the wave function in eq. (fl"2"|) in the 
calculation of the proton electric form- factor in the quark- 
diquark picture of baryon structure. 



3 The proton charge form factor. 

The electric form factor is denned in the standard way [T5] 
as the matrix element of the charge component, Jo(r), 
of the proton electric current between the states of the 
incoming, , and outgoing, p f , electrons in the dispersion 
process, 



q|) =< p/| J (r)|pi >, q = P;-p/- 



(13) 



The explicit evaluation of eq. (flUj) . with Jq{y) — e p |?/>g S t(r)| 2 , 
e p = 1, and plane waves for the electron states amounts 




Q* [CeV 2 ] 

Fig. 3. Comparison of the proton electric form factor as ob- 
tained from the analytic expression in eq. (|15[) (upper curve), 
and the numerical solution of the Bethe-Salpeter equation 
(lower curve ) reported in ref. |16| . In elastic electron-proton 
scattering, Q 2 = — q 2 = —go +<1 2 ! an d go = 0. The experimen- 
tal data set is same as in 1161. 



to the calculation of the following integral, 



|q|) = 



dlr 



(7WH)) 2 sin|q||r| 

MM 



(14) 



In the J* — > J °° limit, the integral calculates exactly 
and is given in terms of q = qd as 



G|(|q|) 



x tan 



b{b 2 + 1) 

|q| 



16|q|b 



q 4 + 4(2& 2 - l)q 2 + 16b 2 (b 2 + 1) 



(15) 



The exact electric form factor of the proton obtained from 
eq. (fTS")) is displayed in Fig. 3 and follows pretty well the 
experimental data. 



4 Effective gluon propagator. 

The proximity of the TQC potential to the QCD quark- 
gluon dynamics is suggestive of the idea to exploit the 
Born approximation and introduce an effective gluon prop- 
agator as a Fourier transform of that very potential. In so 
doing one encounters the integral 



m 



1 2n 



d 3 i 



V(t), 



(16) 



which, unfortunately diverges. As a remedy, we consider 
instead the Fourier transforms of potential matrix ele- 
ments and find integrals that can be taken in closed form. 
Specifically, the Fourier transform of the ground state ma- 
trix element is calculated as 



n (sst) (\q\) = d 



H ' ' ' ^ A f(|r|)^ st (r)r 2 
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Fig. 4. Comparison of the Coulomb propagator (dotted line) 
with the TQC-potential propagator for 6 = and normalized 
ground state wave function (solid line). 
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Fig. 5. Comparison of the Coulomb propagator (dotted line) 
with the TQC-potential propagator for b = 5.85 and unnor- 
malized ground state wave function (solid line). 

, o f n , e~ 2by sin 2 w sin o'w cot w , . 

= Ir*^' <17) 

which expresses in terms of Ei{u) — ^-j— dt, and Ci(u) 
7 + Inu + dt cos *- 1 functions, where 

t,= |a| = M g '^|q| = v ^d. (18) 

In Figs. 4,5 we display 77^ (|q|) for various values of the 
b parameter. The convenience of our consideration is sup- 
ported by the great similarity in the asymptotic behavior 
of the Coulomb propagator and 7l( gst )(|q|). 

5 Concluding remarks. 

In this work we pointed out that the trigonometric Rosen- 
Morse potential captures quite realistically the traits of 



the QCD quark-gluon dynamics. This is reflected by the 
facts that (i) the ratio of (—1/3) of the strengths of the 
linear to Coulomb term following from the Taylor expan- 
sion of the TQC potential is in line with the value of 
(-1/2) - (-1/4) calculated by lattice QCD, (ii) the pre- 
diction of the proton electric form factor in closed form 
follows quite satisfactory data, (iii) the nucleon spectrum 
is well reproduced (and the A one as well, skipped out 
here), (iv) the effective gluon propagator obtained as the 
Fourier transform of the ground state matrix element of 
the TQC potential has reasonable shape and asymptotic 
behavior. The exact single particle basis of the TQC po- 
tential seems quite efficient, indeed, and provides a good 
starting point for more detailed spectroscopic studies. The 
TQC potential can be enriched, for example, by the inclu- 
sion of Oi-Oj interactions to account for the spin dynamics, 
and by screening effects along the line of, say, ref . [3] . The 
effective gluon propagator can be employed in studies of 
non-perturbative QCD phenomena in the spirit of ref. |17j . 
In conclusion, the single particle basis of the TQC poten- 
tial considered here seems to provide an efficient tool for 
quark model calculations of spectroscopic characteristics 
of baryons. 
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